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Abstract 

Low-energy limit of quantum chromodynamics (QCD) is obtained using a mapping theorem recently proved. This 
theorem states that, classically, solutions of a massless quartic scalar field theory are approximate solutions of Yang- 
Mills equations in the limit of the gauge coupling going to infinity. Low-energy QCD is described by a Yukawa theory 
further reducible to a Nambu-Jona-Lasinio model. At the leading order one can compute glue-quark interactions and 
one is able to calculate the properties of the cr and r\ — rf mesons. Finally, it is seen that all the physics of strong 
interactions, both in the infrared and ultraviolet limit, is described by a single constant A arising in the ultraviolet by 
dimensional transmutation and in the infrared as an integration constant. 
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1. Introduction 

Understanding low-energy QCD is a key to uncover 
the structure of light unfiavored mesons. To paraphrase 
Luciano Maiani:"If these are all tetraquark states, where 
are exotic states?"!!]. We will give a first answer to this 
question. 

KLOE-2 measurements[2] seem to rule out quark 
contributions to the structure of lighter mesons. Simi- 
larly, their results support a glue component in 77' meson 
that, having a glue component in the structure, is seen 
to possibly emit a glue state that finally decays in two 
pions. Consistency is only obtained if the emitted glue 
state is a cr meson. 

Structure of cr resonance is hotly debated. Common 
view is that this meson should be a tetraquark state, 
member of a low-lying nonet. This state has not been 
seen yet in lattice computations either quenched or not. 
Recent analysis on yy decay and data from NA48/2[3] 
seem to support the idea that this is a glue state rather 
than a quark composite particle. 

Similarly, f0(980) appears to be a possible glue state 
and can be seen as an excited state of cr. Mass of this 
state should take into account its KK decay. 

All properties of hadronic processes must be consis- 
tent with a single parameter of the theory, A, represent- 
ing a constant arising in ultraviolet limit by dimensional 



transmutation and in the infrared limit as an integration 
constant of the theory. Theory must be working with a 
single A parameter, even if some authors admit a possi- 
ble dependency on the energy scale [4]. Our aim here is 
to fix this constant from a low-energy theoretical analy- 



2. Mapping theorem and low-energy QCD 

In order to manage QCD, it appears essential to find 
a way to reduce this theory to a simpler one. With this 
aim in mind the following theorem has been proved: 

MAPPING THEOREM: An extremum of the ac- 
tion 
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is also an extremum of the SU(N) Yang-Mills La- 
grangian when one properly chooses A" with some com- 
ponents being zero and all others being equal, and A = 
Ng 2 , being g the coupling constant of the Yang-Mills 
held, when only time dependence is retained. In the 
most general case the following mapping holds 
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being if constant, that becomes exact for the Lorenz 
gauge. 

A first proof of this theorem was given in [5] and, 
after a criticism by Terence Tao, a final proof was pre- 
sented in [6] also agreed with Tao||7[]. 
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Applying the above mapping theorem, holding in the 
large coupling limit, QCD generating functional be- 
comes: 

e'/^W-^M+OCl/ViVg). 

At this order ghost field just decouples so we can safely 
ignore it. QCD is so reduced to a Yukawa model by 
the use of the mapping theorem, at the leading order of 
a development in the inverse of the 't Hooft coupling. 
All the parameters of the model are fixed by QCD. In 
the limit A — > oo, the scalar field term takes a Gaussian 
form [8]. This theory is trivial in this limit and the beta 
function goes like /3(A) - A A in four dimensions at lower 
momenta. 

The generating functional for the Yukawa model can 
be cast into a Gaussian form in the strong coupling limit 
g — > oo, knowing the gluon propagator. Indeed, one can 
write 

zmu^ei^-vi^'^'-y^ j n q \.d q wq\x 
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provided y'0 - j ■ rj and the gluon propagator 
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bosons that can be seen as due to Yang-Mills field self- 
interaction. These are the physical states in a strong 
coupling limit. 

We realize that the low-energy limit of QCD can be 
further reduced to a Nambu-Jona-Lasinio model. In the 
gluon propagator we just take the low momenta limit 
producing the contact interaction ilOll 



A(x-y)* — 5\x-y). 
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being 



A first analysis of simple interactions may be accom- 
plished by neglecting quark loops or decays involving 
quarks, just quark-glue vertexes. So, for a first under- 
standing we just consider quark-glue interaction. In this 
case, QCD is exactly integrable producing a non-trivial 
Gaussian generating functional. 

So, the generating functional can be finally integrated 
producing the final result 

■Z[)7,)7,;>]»exp{^(A' 2 -l)/fl' 4 j:c( 4 vv>(j;)A(x-j)j (i (y)}x 
ex P{' fd i xd 4 y 2, n q (x)S [j t ,x-y]ri q (y)} 

The quark propagator, considering that a gradi- 
ent expansion corresponds to a strong coupling 
expansion I llll . is 

S [j>, x - y] = G(t x - t y )6 3 (x -y)x 

(«yroyrl" 11 df Jd i xiA(f-ty-t 11 ,x-y-xi)jf(xi)} 

We have a defined leading order term for a strong cou- 
pling expansion in QCD. Higher order terms contain 
gradients of the fields. QCD at the leading order in a 
strong coupling expansion [g — > 00) appears a confining 
and yet renormalizable theory. 
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We get a first key formula and this is the spectrum of the 
theory, in the strong coupling limit, given by 
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From the mass spectrum we can identify a string tension 

as 



3. cr and 17 - 17' mesons 

We can identify the cr meson as the lowest state in 
Yang-Mills theory. Being massive, this theory shows up 
a mass gap. The mass gap, that is also the mass of the 
lowest glue excitation, is given by 

* = K -(6na s )*A. (8) 
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Width is given by (H (G' L = MS* 
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Presently, the parameter A appears rather arbitrary. Be- 
ing an integration constant, it should be obtained from 
experiment. We just note from this that <Tsu(2)/o~su(3) — 
V2/3 as seen on lattice [9]. We recognize that, at lower 
energies, strong interactions are mediated by a kind of 
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Decay constants for all glue excitations can be 
straightforwardly obtained. Using the mapping theo- 
rem, one has a Fourier series for SU(3) 

^(o,o^ ViSF^r ir=o z," o(-d" ,0-Xlt x^'""'- M + c . c . (io) 



From this series we can easily read the decay constants 
for the glue excitations 



JS " y6na s K(iy ' 1 + e -(2n+l>r • 
Finally, for the <x, setting n — 0, we have 



fa 



cr 2tt 



6na s K(i) 1 + e~" 



(11) 



(12) 



It is interesting to note that, using again the mapping 
theorem, the correlator is given by 

<a;(o, t)A b v (o, o)> = rftf y {tfp, 00(0, o)> = 

oo 

J] B„e- im "'e- w + c.c. + O (l/ V^Vg) 

n=0 

and so, eq.(f5]l is indeed the spectrum of the theory in a 
strong coupling limit. 

77' decay is one of the key processes to understand 
glue role in hadronic physics. This is mainly 77' — > 
■qn + n~ . Measures at DA<BNE proved that 77' has a signif- 
icant glue component [2]. So, we consider this decay as 
a two step process and essentially due to glue emission: 
77' — > j]<r. This will give the following width 
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The opposite process err] — > 77' can be also easily com- 
puted by time reversal symmetry of QCD. Glue produc- 
tion by 77' is a threshold process and so it can be used to 
fix cr mass as 

(13) 



We can use the measured width 11311 to determine f n . 
Taking for r^ the same value observed for the process 
77' — > rpi + if~ and f lT the one computed above, one has 
for 77v at the threshold 



|/,| * 0.019 GeV. 



(14) 



Now, we assume 77 and 77' to mix with an angle 6 ~ -14° 
and also f„ * 0.13 GeV, f ~ -0A5f„ and / 8 * 1.2 f„, 
so one has 



/„ = /o cos + f s sin 9 * -0.019 GeV. 
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while fjf « 0.16 GeV. This computation does not de- 
pend on the value of a s . 



4. QCD and A constant 

Now, we check above scenario against a proper value 
of a s . One has for the running coupling with six 
flavors! 13] 



a s (q , A) 
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and we stop at this order to avoid dependencies from a 
renormalization scheme. Our aim is to fix A in order to 
get a consistent scheme for low-energy QCD. This is the 
same parameter both for high and low-energy physics. 
In our case, a key quantity is string tension. We will 
have 



cr(q 2 ,A) = ^6na s (q 2 , A)A 2 



(17) 



So, the mass of the first glue state is obtained by solving 
the equation 



m. 



= A/0"(777?., A). 



(18) 



2K(i) 
The next excited glue state will be given by 

m\ = —?— Jo-{m 2 ,,h). (19) 
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We can write down the corresponding decay widths as 



^ = -G'^ L (7772,A)777^ 4 Jl - ^ 



and for f0(980) 

2 



(20) 



Tinn = -G' NJL (m 2 , A)mitf Jl f. 

11 V 777T 



(21) 



For the decay constants of glue states we will have 

cr(77^ n ,A) 2n e-("^> 
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On a similar ground we can write, as observed for 
f0(980) from its KK decay, 
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We can compute the ratio rf nK - \gfKK\l\gfnn\ predicted 
to be 2. 59(1. 34) by Mennessier , Narison, and Wang |3J]. 
These authors also agree with a content of f0(980) being 
mostly glue. We have 
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and similarly for decay to Ks. This implies a sizable 
coupling of the cr with K mesons. This computation will 
represent a leading order approximation as we neglect 
mixing effects that should be anyhow present. Finally, 
we will check our computation of A against the corre- 
sponding value of a s as obtained from other sources and 
given by PDG__3]. 
So our estimation is 

A = 0.171 + 0.001 GeV (25) 

that yields the following results for cr mass 

ma- = 0.410 ± 0.007 GeV, (26) 

f0(980) mass 

w /o( 980) = 1 -023 ± 0.002 GeV, (27) 

mass ratio 

m /0 ( 9 80)/«V = 2.49 ±0.05, (28) 

very near the theoretical value 3 but here we are consid- 
ering running coupling and string tension. This may be 
relevant for understanding lattice results, cr width is 



rV/2 = 0.260 ± 0.001 GeV. 



(29) 



It is interesting to note that these values for width and 
mass are very near those of a recent analysis |14|]. Cou- 
pling ratio for f0(980) 



r fnK = 1.42 ±0.02 
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in close agreement with Mennessier , Narison, and 
Wang _3_. cr decay constant is /_. = 0.139 ± 0.008 GeV. 
Finally, we just note that the scenario given through this 
analysis matches rather well that give n in a pioneer- 
ing work of Narison and Veneziano l!l5[ll6ll . With then- 
choice of A, they get a width for cr very near to the 
correct one obtainable from our formulas and a corre- 
sponding decay constant increasing with mass. 

5. Conclusions 

Low-energy limit of QCD can be obtained by an expan- 
sion at very large coupling and remapping Yang-Mills 
field on a quartic massless scalar field theory. This map- 
ping is a proved mathematical theorem. A mass gap 
arises due to strong self-interaction of the fields already 
at classical level. QCD infrared limit is a renormaliz- 
able Yukawa model that, in the proper approximation, 
reduces to a Nambu-Jona-Lasinio model. This is pos- 
sible as we have an analytical formula, in closed form, 



for the gluon propagator in close agreement with lat- 
tice data. A lot of hadronic processes can be described 
by a generating functional obtained in closed form from 
QCD in the above limit but neglecting quark loops and 
quark-quark interaction and having only quark-glue in- 
teractions. From this, one can compute decay constants, 
width and several other experimental meaningful ob- 
servables. cr and f0(980) appear to be glue excitations 
entering into a lot of hadronic processes and describ- 
ing the true physical states of glue for low-energy QCD. 
It is worthwhile to note that decay constants for glue 
excitations are strongly damped for higher excitations 
making observable just the first few states. r( appears 
mostly to emit a glue state that, due to a threshold con- 
dition, should properly fix the value of the mass of cr 
resonance. In agreement with KLOE-2 evidences, this 
gives a sound explanation of the main decay of this me- 
son resembling the case of QED of an excited atom de- 
caying into an atom plus the force carrier, i.e. a photon. 
This decay can also be used to determine the decay con- 
stant of r]. All observables are properly fixed with a sin- 
gle constant A = 171 ± 1 MeV that recovers the correct 
value of the measured running coupling at Mz- 
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